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Abstract: Various performance indices are used for the design of serial 
manipulators. One method of optimization relies on the condition num- 
ber of the Jacobian matrix. The minimization of the condition number 
leads, under certain conditions, to isotropic configurations, for which the 
roundoff-error amplification is lowest. In this paper, the isotropy condi- 
tions, introduced elsewhere, are the motivation behind the introduction 
of isotropic sets of points. By connecting together these points, we define 
families of isotropic manipulators. This paper is devoted to planar manipu- 
lators, the concepts being currently extended to their spatial counterparts. 
Furthermore, only manipulators with revolute joints are considered here. 

1. Introduction 

Various performance indices have been devised to assess the kinetostatic 
performance of serial manipulators. The literature on performance indices 
is extremely rich to fit in the limits of this paper, the interested reader 
being invited to look at it in the rather recent references cited here. A di- 
mensionless quality index was recently introduced by Lee, Duffy, and Hunt 
(1998) based on the ratio of the Jacobian determinant to its maximum ab- 
solute value, as applicable to parallel manipulators. This index does not 
take into account the location of the operation point in the end-effector, for 
the Jacobian determinant is independent of this location. The proof of the 
foregoing fact is available in (Angeles, 1997), as pertaining to serial manip- 
ulators, its extension to their parallel counterparts being straightforward. 
The condition number of a given matrix, on the other hand is well known 



to provide a measure of invertibility of the matrix (Golub and Van Loan, 
1989). It is thus natural that this concept found its way in this context. 
Indeed, the condition number of the Jacobian matrix was proposed by Sal- 
isbury and Craig (1982) as a figure of merit to minimize when designing 
manipulators for maximum accuracy. In fact, the condition number gives, 
for a square matrix, a measure of the relative roundoff-error amplification of 
the computed results (Golub and Van Loan, 1989) with respect to the data 
roundoff error. As is well known, however, the dimensional inhomogeneity 
of the entries of the Jacobian matrix prevents the straightforward applica- 
tion of the condition number as a measure of Jacobian invertibility. The 
characteristic length was introduced in (Angeles and Lopez-Cajun, 1992) 
to cope with the above-mentioned inhomogeneity. Apparently, neverthe- 
less, this concept has found strong opposition within some circles, mainly 
because of the lack of a direct geometric interpretation of the concept. It 
is the aim of this paper to shed more light in this debate, by resorting to 
the concept of isotropic sets of points. Briefly stated, the application of 
isotropic sets of points to the design of manipulator architectures relies on 
the concept of distance in the space of m x n matrices, which is based, in 
turn, on the Frobenius norm of matrices. With the purpose of rendering 
the Jacobian matrix dimensionally homogeneous, moreover, we introduce 
the concept of posture-dependent conditioning length. Thus, given an arbi- 
trary serial manipulator in an arbitrary posture, it is possible to define a 
unique length that renders this matrix dimensionally homogeneous and of 
minimum distance to isotropy. The characteristic length of the manipula- 
tor is then defined as the conditioning length corresponding to the posture 
that renders the above-mentioned distance a minimum over all possible 
manipulator postures. 

It is noteworthy that isotropy comprising symmetry at its core, manip- 
ulators with only revolute joints are considered here. It should be apparent 
that mixing actuated revolutes with actuated prismatic joints would destroy 
symmetry, and hence, isotropy. 

2. Algebraic Background 

When comparing two dimensionless m x n matrices A and B, we can 
define the distance d(A, B) between them as the Frobenius norm of their 
difference, namely, , 



d(A, B) = ||A - B|| i.e., d(A, B) = W-tr[(A - B)(A - B) T ] (1) 



An mx n isotropic matrix, with m < n, is one with a singular value a > 
of multiplicity m, and hence, if the m x n matrix C is isotropic, then 




CC T = a 2 l 



(2) 



where 1 is the rax m identity matrix. Note that the generalized inverse of 
C can be computed without roundoff-error, for it is proportional to C T , 
namely, 



(CC 7 )" 1 ^ = -ic 
a 1 



(3) 



Furthermore, the condition number k(A) of a square matrix A is defined 
as (Golub and Van Loan, 1989) 



\A\\\\A- 



(4) 



where any norm can be used. For purposes of the paper, we shall use the 
Frobenius norm for matrices and the Euclidean norm for vectors. Hence- 
forth we assume, moreover, a planar n-revolute manipulator, as depicted 
in Fig. 1, with Jacobian matrix J given by (Angeles, 1997) 
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(5) 




where r, is the vector directed from the center of the ith revolute to the 
operation point P of the end-effector, and matrix E represents a counter- 
clockwise rotation of 90° . 

It will prove convenient to partition J into 
a 1 x n block A and a 2 x n block B, defined 
as A = [1 1 • • • 1] and B = [Eri Er 2 • • • Er n ]. 
Therefore, while the entries of A are dimension- 
less, those of B have units of length. Thus, the 
sole singular value of A, i.e., the nonnegative 
square root of the scalar of AA T , is \fn, and 
hence, dimensionless, and pertains to the map- 
ping from joint-rates into end-effector angular 
velocity. The singular values of B, which are the 
nonnegative square roots of the eigenvalues of 

BB T , have units of length, and account for the mapping from joint-rates 
into operation-point velocity. It is thus apparent that the singular values 
of J have different dimensions and hence, it is impossible to compute re(J) 
as in eq.(4), for the norm of J cannot be defined. The normalization of the 
Jacobian for purposes of rendering it dimensionless has been judged to be 
dependent on the normalizing length (Paden, and Sastry, 1988; Li, 1990). 
As a means to avoid the arbitrariness of the choice of that normalizing 
length, the characteristic length L was introduced in (Ranjbaran, Angeles, 
Gonzalez-Palacios, and Patel, 1995). We shall resort to this concept, while 
shedding more light on it, in discussing manipulator architectures. 



Figure 1. Planar n-revolute 
manipulator 



3. Isotropic Sets of Points 



Consider the set S = { Pk }" of n points in the plane, of position vectors 
{ pk }", and centroid C, of position vector c, i.e., 



1 n 

= -Y]pfc 



" 1 



(6) 



The summation appearing in the right-hand side of the above expres- 
sion is known as the first moment of S with respect to the origin O from 
which the position vectors stem. The second moment of S with respect to 
C is defined as a tensor M, namely, 

n 

i 

It is now apparent that the root-mean square value of the distances 
{dk}™ of S, d Tms , to the centroid is directly related to the trace of M, 
namely, 



Otitis — 



lX>*-c) T (P*-c)^tr(M) 



^ n i 



(8) 



Further, the moment of inertia I of S with respect to the centroid is 
defined as that of a set of unit masses located at the points of S, i.e., 

n 

I = J} IIP* " c|| 2 l - (p fc - c)(p fc - cf } (9a) 
l 

in which 1 is the 3x3 identity matrix. Hence, in light of definitions (7) and 
^' I = tr(M)l - M (9b) 

We shall refer to I as the geometric moment of inertia of S about its 
centroid. It is now apparent that I is composed of two parts, an isotropic 
matrix of norm tr(M) and the second moment of S with the sign reversed. 
Moreover, the moment of inertia I can be expressed in a form that is more 
explicitly dependent upon the set {p/c — c}™, if we recall the concept of 
cross-product matrix (Angeles, 1997): For any three-dimensional vector v, 
we define the cross-product matrix of (p& — c), or of any other three- 
dimensional vector for that matter, as 

P* - a " Pt ~ V C) X V| (10a) 



Further, we recall the identity (Angeles, 1997) 

p| = -Hp* - c|| 2 i + ( Pfc - c)( P * - c) T (iob) 

It is now apparent that the moment of inertia of S takes the simple 
form n 

i = -E p * (ii) 
i 

We thus have 

Definition 1 (Isotropic Set) The set S is said to be isotropic if its second- 
moment tensor with respect to its centroid is isotropic. 

As a consequence, we have 

Lemma 1 The geometric moment of inertia of an isotropic set of points 
is isotropic. Conversely, an isotropic geometric moment of inertia pertains 
necessarily to an isotropic set of points. 

We describe below some properties of isotropic sets of points that will help 
us better visualize the results that follow. 



3.1. ISOTROPY-PRESERVING OPERATIONS ON SETS OF POINTS 

Consider two isotropic sets of points in the plane, Si = { Pk }" and S2 = 
{ Pk }n+T- If the centroid C of the position vector c of Si coincides with 
that of 52, i.e. if, 

^ n ^ n+m 

c = -Vpfc = — Vpfc (12) 
n ^ m ^ 

1 ra+l 

then, the set S = Si U £2 is isotropic. Hence, 

Property 1 The union of two isotropic sets of points sharing the same 
centroid is also isotropic. 

Furthermore, as the reader can visualize, we state below one more 
operation on sets of points, namely, a rigid-body rotation, that preserves 
isotropy: 

Property 2 The rotation of an isotropic set of points as a rigid body with 
respect to its centroid is also isotropic. 



3.2. TRIVIAL ISOTROPIC SETS OF POINTS 

An isotropic set of points can be defined by the union, rotation, or a com- 
bination of both, of isotropic sets. The simplest set of isotropic points is 
the set of vertices of a regular polygon. We thus have 



Definition 2 (Trivial isotropic set) A set of n points S is called trivial 
if it is the set of vertices of a regular polygon with n vertices. 

Lemma 2 A trivial isotropic set S remains isotropic under every reflection 
about an axis passing through the centroid C . 

4. An Outline of Kinematic Chains 

The connection between sets of points and planar manipulators of the serial 
type is the concept of simple kinematic chain. For completeness, we recall 
here some basic definitions pertaining to this concept. 

4.1. SIMPLE KINEMATIC CHAINS 

The kinematics of manipulators is based on the concept of kinematic chain. 
A kinematic chain is a set of rigid bodies, also called links, coupled by 
kinematic pairs. In the case of planar chains, two lower kinematic pairs are 
possible, the revolute, allowing pure rotation of the two coupled links, and 
the prismatic pair, allowing a pure relative translation, along one direction, 
of the same links. For the purposes of this paper, we study only revolute 
pairs, but prismatic pairs are also common in manipulators. 

Definition 3 (Simple kinematic chain) A kinematic chain is said to be 
simple if each and every one of its links is coupled to at most two other links. 

A simple kinematic chain can be open or closed; in studying serial 
manipulators we are interested in the former. In such a chain, we distinguish 
exactly two links, the terminal ones, coupled to only one other link. These 
links are thus said to be simple, all other links being binary. In the context 
of manipulator kinematics, one terminal link is arbitrarily designated as 
fixed, the other terminal link being the end-effector (EE), which is the one 
executing the task at hand. The task is defined, in turn, as a sequence of 
poses — positions and orientations — of the EE, the position being given at 
a specific point P of the EE that we term the operation point. 

4.2. ISOTROPIC KINEMATIC CHAINS 

To every set S of n points it is possible to associate a number of kinematic 
chains. To do this, we number the points from 1 to n, thereby defining n — 1 
links, the iih link carrying joints i and i + Links are thus correspondingly 
numbered from 1 to n, the nth link, or EE, carrying joint n on its proximal 
(to the base) end and the operation point P on its distal end. Furthermore, 
we define an additional link, the base, which is numbered as 0. 

It is now apparent that, since we can number a given set S of n points 
in n\ possible ways, we can associate n\ kinematic chains to the above set S 



of n points. Clearly, these chains are, in general, different, for the lengths of 
their links are different as well. Nevertheless, some pairs of identical chains 
in the foregoing set are possible. 

Definition 4 (Isotropic kinematic chain) Let set S of n points be iso- 
tropic, and the operation point P be defined as the centroid of S. Any kine- 
matic chain stemming from S is said to be isotropic. 



5. The Posture-Dependent Conditioning Length of Planar n-Re- 
volute Manipulators 

Under the assumption that the manipulator finds itself at a posture V that 
is given by its set of joint angles, { Ok }™ , we start by dividing the last two 
rows of the Jacobian by a length l-p, as yet to be determined. This length 
will be found so as to minimize the distance of the normalized Jacobian to 
a corresponding isotropic matrix K, subscript V reminding us that, as the 
manipulator changes its posture, so does the length l-p. This length will be 
termed the conditioning length of the manipulator at V. 



5.1. A DIMENSIONALLY-HOMOGENEOUS JACOBIAN MATRIX 

In order to distinguish the original Jacobian matrix from its dimensionally- 
homogeneous counterpart, we shall denote the latter by J, i.e., 

1 1 ••• 1 

(1/h) Eri (l/l v ) Er 2 • • • (l/l v ) Er„ 

Now the conditioning length will be defined via the minimization of 
the distance of the dimensionally-homogeneous Jacobian matrix J of an n- 
revolute manipulator to an isotropic 3xn model matrix K whose entries are 
dimensionless and has the same gestalt as any 3 x n Jacobian matrix. To this 
end, we define an isotropic set K. = { Ki }™ of n points in a dimensionless 
plane, of position vectors { kj }™, which thus yields the dimensionless matrix 



K 



1 1 
Eki Ek 2 



1 

Ek„ 



(13) 



Further, we compute the product KK r : 

EikfE^ 



KK 



n 



E"Ek t EiEk,kfE r J 



Upon expansion of the summations occurring in the above matrix, we 



have 



]TkfE r = (£Ek 4 f = E(J>f and ^Ek.kfE^ = E(J>kf)E^ 



Now, by virtue of the assumed isotropy of /C, the terms in parentheses 
in the foregoing expressions become 



= and ^kjkf = /c 2 l 2x2 
i l 

where the factor k 2 is as yet to be determined and 1-2x2 denotes the 2x2 



identity matrix. Hence, the product KK T takes the form 



KK 



(14) 



n T 
fc 2 l 2x2 

Now, in order to determine k 2 , we recall that matrix K is isotropic, and 
hence that the product KK T has a triple eigenvalue. It is now apparent 
that the triple eigenvalue of the said product must be n, which means that 

n 

k 2 = n and ^k^kf = (n)l 2X 2 (15) 
l 

5.2. COMPUTATION OF THE CONDITIONING LENGTH 

We can now formulate a least-square problem aimed at finding the condi- 
tioning length l-p that renders the distance from J to K a minimum. The 
task will be eased if we work rather with the reciprocal of lp, A = I /l-p, 
and hence, 

z = --tr[(J -K)(J -K) T ] -► min (16) 
Z ft A 

Upon simplification, and recalling that tr(KJ ) = tr(JK T ), 

z = --tr(JJ T -2KJ T + KK T ) (17) 
2 n 

It is noteworthy that the above minimization problem is (a) quadratic 
in A, for J is linear in A and (b) unconstrained, which means that the 
problem accepts a unique solution. This solution can be found, additionally, 
in closed form. Indeed, the optimum value of A is readily obtained upon 
setting up the normality condition of the above problem, namely, 

3z 1 (d{3~f)\ 1 („df\ n 

M = fc te (T)-; b ( K -5r) =0 (18) 

where we have used the linearity property of the trace and the derivative 
operators. The normality condition then reduces to 



Now, if we notice that ||rj || is the distance dj from the operation point P 
to the center of the jth revolute, the first summation of the above equation 
yields nd^ ms , with d rms denoting the root-mean-square value of the set of 
distances { dj }", and hence, 



A 



nd^ ms 



h 



n( ^rms 



(19) 



Definition 5 (Characteristic length) The conditioning length of a ma- 
nipulator with Jacobian matrix capable of attaining isotropic values, as per- 
taining to the isotropic Jacobian, is defined as the characteristic length of 
the manipulator. 

An example is included below from which the reader can realize that a 
reordering of the columns of K preserves the set of manipulators resulting 
thereof, and hence, K is not unique. Lack of space prevents us from includ- 
ing an example of a nonisotropic manipulator. The interested reader can 
find such an example in (Chablat and Angeles, 1999.) 

5.3. EXAMPLE: A FOUR-DOF REDUNDANT PLANAR MANIPULATOR 

An isotropic set K, of four points, {Ki}f, is defined in a nondimensional 
plane, with the position vectors kj given below: 



,k 3 



1 



(20a) 



which thus lead to 



K = 



1 1 
1 1 
-1 1 



KK T = 



4 
4 
4 



= a 2 l 



(20b) 



and is apparently isotropic, with a triple singular value of k = 2. We thus 
have 4! = 24 isotropic kinematic chains for a four-dof planar manipulator, 
but we represent only six in Fig. 2 because the choice of the first point is 
immaterial, since this choice amounts to a rotation of the overall manip- 
ulator as a rigid body. The conditioning length l-p is equal to 1/2 for the 
three manipulators, with / denoting any nonzero length. Moreover, the link 
lengths are defined, respectively, for three distinct manipulators of Fig. 2, 
form left to right: a\ = I, a 2 = I, a% = I, a^ = \/2Z/2; a\ = I, a 2 = \/2Z, 03 = 
I, a 4 = V21/2; a x = \/2l, a 2 = /, a 3 = V21, a 4 = V21/2. 

An extensive discussion of isotropic sets of points and the optimum 
design of manipulators is available in (Chablat and Angeles, 1999.) 




Figure 2. Six isotropic postures for the same isotropic set 
6. Conclusions 

Isotropic manipulators are defined in this paper by resorting to the con- 
cept of isotropic sets of points. This concept allows us to define families 
of isotropic redundant manipulators. The characteristic length can thus be 
defined as the conditioning length l-p that pertains to the isotropic postures. 
The paper focuses on planar manipulators, but the underlying concepts are 
currently extended to their three-dimensional counterparts. 
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